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Abstract 

Nucleon-nucleon, axion bremsstrahlung is the primary mechanism for axion emission from 
the nascent neutron star associated with SN 1987 A, and the rate for this process has been 
calculated in the one pion exchange approximation (OPE). The axion mass limit which 
follows from SN 1987A, m a < 10 _3 eV, is the most stringent astrophysical bound, and has 
received much scrutiny. It has been suggested that by using OPE to calculate the cross 
section for the analog process, pp — » pp + 7r°, and comparing the result to experimental 
data one can test the validity of this approximation, and further, that such a comparison 
indicates that OPE leads to a value for this cross section which is a factor of 30-40 too 
large. If true, this would suggest that the axion mass limit should be revised upward by a 
factor of ~ 6. We have carefully evaluated the cross section for pp — ► pp + using OPE 
and find excellent agreement (to better than a factor of 2) with the experimental data. 
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I. Introduction 

There has been a great deal of interest in axion emission from SN 1987A. 1-5 And 
indeed, consideration of the effect of axion emission on the neutrino burst observed by the 
KII 6 and IMB 7 detectors seems to provide strong evidence against the existence of an axion 
with mass in the range 10 -3 — 2 eV; for the DFS type axion this improved the existing 
astrophysical bound by a factor of ~ 10, while for the hadronic type axion, the improvement 
was more than a factor of 10 3 . 8 If the axion exists, then the dominant emission process 
from SN 1987A should have been nucleon-nucleon, axion bremsstrahlung. The matrix 
element squared for this process has been computed in the OPE approximation; 2 the 4 
direct and 4 exchange diagrams are shown in Fig. 1. Given that the pion-nucleon coupling 
is of order unity one might question the accuracy of such an approximation — of course, 
one should remember that the pion-nucleon coupling is derived by a comparison between 
an OPE treatment of pion-nucleon scattering and experimental data. 9 In addition, since 
the densities that existed at the core of the nascent neutron star associated with SN 1987A 
shortly after collapse were ~ 2p nuc i ea r ~ 8 x 10 14 g cm -3 , one should also worry about 
collective nuclear effects. Here, we will restrict our discussion to the validity of the OPE 
approximation itself. 

Choi, et al 5 recently suggested a clever way of checking the accuracy of the OPE approx- 
imation in computing nucleon-nucleon, axion bremsstrahlung. Their idea is to compute 
the cross section for the analog process, nucleon-nucleon, pion bremsstrahlung using OPE 
and to compare to the body of existing experimental data. Both the axion and the pion 
are pseudo-Nambu-Goldstone bosons, and as such couple derivatively: 

£i„t = - + i^d^aN^N + + ... (1) 

Am Am 

where N is the nucleon field, 7r° is the neutral pion field, a is the axion field, m ~ 0.94 
GeV is the nucleon mass, g a N ~ m/(fa/N a ) is the axion-nucleon coupling, and in the 
OPE approximation the pion-nucleon coupling A ~ 2 m/m„ ~ m/f n , where m, ~ 135 
MeV is the neutral pion mass and f„ ~ 95 MeV is the pion decay constant. We should 
mention that the derivative coupling of the pion to the nucleon follows directly from the 
chiral Lagrangian, a low-energy, effective Lagrangian which describes the interactions of 
pions and nucleons. The similarity of the pion-nucleon and axion-nucleon couplings, as 
well as the subtle issue of when the pseudo- vector coupling can be expressed as a pseudo- 
scalar coupling, are addressed in Ref. 5. Here, we will use the unambiguously correct 
pseudo- vector coupling for both the axion and pion. 

The axion mass, m a , and PQ symmetry breaking scale, /„, are related by 

m a = 0.62 eV [10 7 GeV/(/ a /JV a )] 

where N a is the color anomaly of the PQ symmetry. For details about the axion and 
its couplings to matter, see Refs. 10. Because of the similarity of the pion and axion 
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couplings we see that by substituting A — ► g a n and m a —* m *■, the matrix element for 
NN —* NN + 7r° can be obtained from that for NN — ► NN + a. 

The matrix element squared for nucleon-nucleon, axion bremsstrahlung has been cal- 
culated in Ref. 2; for the process pp — > pp + a it is 


E ML.. 

SPIN 


256 ffq p ra 2 
3 m\ 


( 3 -/ 9 ) 


( 2 ) 


where |-M| 2 xion has been summed over ioifc, initial and final proton spins and averaged over 
the directions of the axion and the 3-momentum exchanged. The quantity 0 is related to 
the average of the cosine squared of the angle between the direction of the momentum 
transfer in the direct and exchange diagrams: for degenerate matter 0 = 0; and for non- 
degenerate matter 0 (see Ref. 2). For our purposes here we take 0 ~ 1. In addition, 
several approximations were made in calculating \Nl 1 axion- the nucleons were assumed to 
be non-relativistic, the axion mass was taken to be zero, and 3 mT was assumed to be much 
larger than ml (here T = temperature) so that the pion mass in the pion propagator can 
be neglected — all good approximations in the core of the neutron star associated with SN 
1987A. However, we must keep them in mind to understand the realm of validity for using 
| Ad | axion compute the corresponding matrix element squared for pp — + pp + 7r°. 


II. Naive OPE 


Once \M\ 2 n is at hand it is straightforward to obtain the cross section for pp 
(Ref. 11): 


PP + 7T° 


dc = \ S \M\l(2C 12 )-\2Tr) 4 6 4 (p 1 +p 2 -p 3 -p 4 -a)<m z <m A dn a (3) 

SPIN 

where!, 2 denote the incoming protons, 3,4 denote the outgoing protons, a denotes the 
pion, dUi = d 3 pi/2Ei(27r) 3 , S = 1/2 is the usual symmetry factor for identical particles in 
the final state, the factor of 1/4 is inserted to average over initial proton spins, and the 
kinematical factor (Moeller flux factor) 

C 12 = s 1/2 (s — 4 m 2 ) 1 / 2 = s(l — 4 m 2 /s) 1 / 2 

The quantity s is the center-of-mass (CM) energy squared, which is related to the KE of 
the incoming proton in the lab (= Tl) by 

s = 4m 2 -f 2mTi 

The momenta of the incoming protons in the CM frame is 


\ P \\ 2 = \p2 \ 2 = tuTl/2 
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The threshold for pion production is 

Tl (threshold = 2 m*. + m 2 /2m = 280 MeV ~ 2m*. 

Because the axion mass was taken to be zero in computing |.M| 2 xion ^ i s n °t possi- 
ble to obtain \M\\ by the substitution: g ap — ► A. However, one can, as a first, naive 
approximation, do so. Of course, by so doing one is implicitly assuming that the pion is 
ultra-relativistic, i.e., Tl Tj,| t hreshold, so that the approximation m*. = 0 is a good one. 
In the next Section we will compute |.M| 2 without recourse to this assumption; for the 
moment we will use this naive approximation to compute a(pp pp + ir°). From the last 
two formulae we see that for the approximations made in obtaining \A4 | 2 from |Af| 2 xion 
to be valid (relativistic pion, non- relativistic nucleons) Tl must satisfy: 

2 m*. < T l < 2m or 0.3 GeV < T L < 2.0 GeV 

To actually compute the cross section it is most convenient to use the Dalitz represen- 
tation for the 3-body phase space factor: 11 


(27r) 4 ^ 4 (pi + p 2 — pz — Pa — a)dU 3 dn.4dU a = (16s) 1 (27 t) 3 drn\ i drn\ a (4) 

where we have taken advantage of the fact that the matrix element squared is constant to 
perform most of the integrals, and the invariant mass variables are 

m 34 = ( P3 + Pif rn\ a = (p 3 + a) 2 


It is then straightforward to obtain cr(pp — > pp + i r°): 

^pp— . pp + 0 = _l_ £ 


SPIN 


r(y/s — rn x ) 2 

1= I (x — 4m 2 ) 1/,2 [(s — x — m 2 ) 2 — 4xm 2 ] 1//2 -^= 

J ir, 


(5a) 

( 56 ) 


jl/2 dx 

i i m 2 yfl 

where x = m 2 4 . Using the matrix element squared which is obtained from |Af | 2 x ; on by 
substituting g ap — * m/f n we find 


a = 2.46 x 10 25 cm 2 


s 2 (l — Am 2 Is) 1 ! 2 


( 6 ) 


Note, that although approximations have been made in evaluating \M | 2 , the phase space 
integrals have been evaluated exactly. 

In the limit that the 2 outgoing protons and pion are non-relativistic the integral I can 
be evaluated in closed form (see Byckling and Kajantie 11 ), 


I = 


7T (mjr/m) 1 / 2 

y/2(l + m*-/2m) 3 / 2 


s 2 (l — 2m/ -y/s — m*-/ -\/s) 2 


(7) 


4 


* 
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While the assumption that the outgoing protons are non-relativistic is probably reasonable, 
the assumption that the outgoing pion is non-relativistic is not a good one; moreover, it 
is contrary to the approximation made in computing |Af| 2 . However, over the lab energy 
range, 0.3 GeV < Tl <1.0 GeV, this approximate formula for I is suprisingly accurate, to 
about 5%. Using this expression for I, we obtain 


= 1 (m^/m) 1 / 2 ^ 2 (1 - 2m/ - rn^/sf 

2048^2(1 + m,/2m) 3 / 2 ^ 1 '* (1 - 4m 2 / 5 ) 1 /2 ^ ; 

(&b) 

a form similar to that of Choi, et al 5 . 

There exists a wealth of experimental data for the process pp — ► pp -f tt° at energies 
Tl It 300 MeV. The data is summarized in Refs. 12 and 13, and a convenient analytical fit 
is given in Ref. 12. Both the experimental data and the analytical fit are shown in Fig. 2. 
The analytical fit is given by the following expression: 


where 


_ = * / > g m 2 r 2 (g/q 0 ) 3 

2|p a |2 (< M >2 —171q) 2 + m2f2 


-2 ( -\ I s - ( m ~ < M >) 2 ][s - (m+ < M >)2] 

p r {S) - — 

2 ( M 2 x _ [< M > 2 -(m - m„)2][< M > 2 -(m + m^) 2 ] 
^ V ; 4 < M > 2 

9o = firnl) 

pi — (m + mo) 2 / 4 - m 2 

< M >= Mo + (tan " 1 Z + - tan " 1 Z .)" 1 ^ In 

4 1 + Zf. 

Z+ = 2(y/s — m — Mo) /T 0 Z_ = 2(m + m n — Mo) /To 

a = 3.772 0 = 1.262 

Mo = 1.22 GeV m 0 = 1.188 GeV 
r 0 = 0.12 GeV r = 0.099 GeV 


( 9 ) 


The expression we obtained for the cross section by substituting g ap — ► A, Eq.(6), along 
with the experimental data is shown in Fig. 3. Recall that the range of validity for our 
calculation above was: 0.3 GeV •< Tl -C 2 GeV. In that region the experimental data 
agree very well with our calculation, to better than a factor of 2. As one moves outside 
this range our calculation begins to overestimate the cross section significantly. At small 
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Ti our calculation is not valid as the pion is not highly-relativistic: Because the pion cou- 
ples derivatively \M\\ should vanish at threshold, and thus our approximation to \M\\, 
which is independent of Tl, should overestimate |.M|* . (As we shall see shortly, \M.\\ does 
not actually vanish at threshold, rather is suppressed by a factor of m*/m.) Likewise, 
at large Tl, where the nucleons are becoming relativistic, we would expect our approxi- 
mation to overestimate a because the momentum dependence of nucleon propagators has 
been neglected. Thus, the deviation of the OPE treatment of the cross section from the 
experimental data, at both small and large Tl, is as expected. 

In order to estimate the effect of pion production threshold on the OPE calculation of 
a(pp —* pp- 1- tt°) we have calculated |A4|* at threshold (Tl = 2m*) with m a = m* ^ 0. 
(Of course, the threshold dependence of phase space has already been properly taken into 
account.) We find 


E 

SPIN 


m 


l^l. = 64A 2 — , 


mi 


or a factor of 3 m*/8m ~ 1/20 smaller than the value obtained taking m a = 0. 14 Moreover, 
slightly above threshold (2E > 300 MeV), we find that the square of any given diagram is 
modified from the m a = 0 result by a factor 


1 - m 2 m 2 /(p 3 ■ af oc \Scm\ 2 °c (Tt/2m, - 1) 


( 10 ) 


The effect of correcting the naive OPE cross section by this approximate threshold factor 
for the matrix element is also shown in Fig. 3: The agreement between the OPE result and 
the experimental data improves dramatically for small values of Tl. Of course, to properly 
take into account the fact that m* 0 one should compute the complete matrix element 
squared with m* ^ 0 which we will do in the next section. But first let us compare to the 
results of Choi, et al 5 . 

Choi, et al 5 basically carried out the same calculation as we did to obtain the ‘naive 
OPE’ prediction: that is, they used the matrix element squared obtained from axion 
bremsstrahlung and substituted g ap — ► A. They then compared this result to similar 
experimental data at energies Tl = 400, 500 MeV (although there seems to be a slight 
discrepancy between the data they use at 400 MeV, and that in Refs. 12 and 13). They 
found that the OPE approximation overestimated the experimental cross section by factors 
of 30-40. Our disagreements with the calculation of Choi, et al 5 are manifold. Firstly, the 
value of S X^SPIN 1-^1 * that they used is larger than ours by a factor of 12(/*/ra*) 2 ~ 6.6. 
Of this, a factor of 4(/*/m*) 2 ~ 2.0 traces to their using A = 2 m/m* for the pion nucleon 
coupling vs. our using A = m/ /*; the remainder, a factor of 3.3, apparently has to do with 
spin averaging and the symmetry factor S. With regard to the value of A; one can make 
arguments for either choice, and in any case the discrepancy due to this is only a factor of 
2. (The issue of the pion nucleon coupling is addressed further in the Appendix.) Secondly, 
they have apparently used the non-relativistic approximation to evaluate I. However, the 
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numerical factor in their expression for I differs from ours by a factor of 

1 + m rr /2m) 3 / 2 ~ 2.63 

7T 

Taken together, these factors account for an overall factor of 6.58 • 2.63 ~ 17.3, which is 
the entire discrepancy between our expressions for cr(pp — ► pp -f tt°) and their expression. 
Further, they have only compared to the experimental data at two energies, Tl = 400, 500 
MeV, energies which are not too far from the threshold for pion production, where the 
approximation made in using |Af| 2 xion to obtain \M | 2 (i.e., relativistic pion) is breaking 
down. As mentioned earlier, near threshold one expects this simple technique to overesti- 
mate cr(pp — + pp + 7T°). And as we will see shortly, when the exact matrix element is used 
the discrepancy at threshold disappears. 

III. Exact OPE 


We have calculated the full matrix element squared for the process pp — » pp -f tc° in 
the OPE approximation, making no other approximations: that is, all particles treated 
fully relativistically and the mass of the bremsstrahlung pion no longer taken to be zero 
(the square of the sum of the diagrams in Fig. 1 with the axion replaced by a neutral 
pion). The resulting expression for |Af| 2 is considerably more complicated than Eq.(2), 
and is given in the Appendix. The OPE cross section for pp — » pp + n° was then computed 
by Monte Carlo integration of the matrix element squared over phase space. Three-body 
phase space is characterized by 4 independent parameters which we take to be: p 3 , p 4 , the 
3-momenta of the outgoing protons in the CM frame; and 73, 74, the cosine of the angles 
between the outgoing protons and incoming protons, again in the CM frame. The cross 
section is then given by: 


a = 


1 m 2 (m^/m) 1 / 2 (1 — 2 m/y/s — m„/y/s ) 2 f \Ai\ 2 dP 

512y^7r 2 m*f 2 (1 + m n /2m) 3 / 2 (1 — 4m 2 /s) 1 / 2 f dP 

dP = £e!1 cos 4,(1 - T j)U(i - T j )■/. * >*»*»*»■ 


( 11 ) 


I 


7 r 2 (m n lrrij l l 2 , n 

iP = ^ W2m)V’ (1 ~ 2m/v ^ " m ' /v ^ 


sin <f> = 


7374 — 734 


(1 - 73 2 ) 1/2 (1 -74 2 ) 1/2 


s + 2m 2 - m 2 + 2E 3 Ei - 2 v /s(£ 3 + E4) 

734 = 

2P3P4 

where <f> is the azimuthal angle of outgoing proton 4 (around the beam axis), 734 is the 
cosine of the angle between the two outgoing protons, and M = (2 m 2 / f n m\)M.. In 
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addition, the total volume of the 3-body phase space (= /dP) has been approximated by 
Eq.(7) (which we recall is accurate to about 5%). The full matrix element squared \M\ 2 
is given in the Appendix. While the 4 integration variables are independent, there are 
complicated kinematical constraints amongst them, i.e., the boundaries of the region of 
integration are not easily expressed in closed form. We have implemented the constraints 
in the following way: the 4 variables p 3 , p 4 , 73, 74 are picked at random, and then checked 
to see if they satisfy the constraints; if they do, then the matrix element is evaluated, and 
the contribution to the integral computed. For each Monte Carlo integration we selected 
10 4 sets of the 4 variables, of which about 3 x 10 3 satisfied the kinematical constraints; 
thus one expects the integration to have an accuracy of a few percent. Our results are 
shown in Fig. 4, along with the experimental data. 

The agreement between the ‘exact’ OPE result and the data is quite impressive. The 
only significant deviation is in the region Tjr, ~ 500 — 900 MeV where the OPE result 
underestimates cr by a factor of up to 3. The reason for the underestimation is simple to 
understand: The threshold for A(1232) production is 630 MeV, and in the region of 500- 
800 MeV the contribution of the 3,3 resonance (pp — ► Ap — > pp-f- 7r°) should be important, 
and is not accounted for in the OPE approximation. (In addition, there is the 1V , (1440) 1,1 
resonance whose production threshold is 1130 MeV.) In fact, with a little imagination one 
can see the Breit-Wigner shape in the deviation between the OPE results and the data. 

As a check, we have compared the exact OPE result above, with the pion mass in 
the matrix element squared set to zero (but not in the phase space), with the previous 
naive OPE result, i.e., Eq.(8b). The comparison is shown in Fig. 5. The previous results 
are reproduced with good accuracy except near threshold. In fact, this can be easily 
understood: In using Eq.(2) for the matrix element squared we have set /? = 1; exactly 
at threshold the 3-momentum transfer in the direct and exchange diagrams are equal and 
—y 3. Had we used f3 = 3 in evaluating Eq.(2) the matrix element squared would 
vanish.' Of course, when the matrix element squared is evaluated exactly, the angular 
average corresponding to /3 is automatically calculated. Further, as an additional check, 
we have compared the exact result (with m n ^ 0) to the naive result at pion production 
threshold. As discussed earlier, in the naive approximation the matrix element squared is 
overestimated by a factor of ~ 20 at threshold, which is what we find. 

IV. Discussion/Summary 

In the supernova, where the nucleons have thermal distributions characterized by tem- 
peratures of order 20-80 MeV, the thermally averaged CM energy is given by 

< s >=< ( pi + P2) 2 >— < 4m 2 + |pi| 2 + |p2| 2 — 2pi • P2 >~ 4m 2 + 6mT 

so that the average value of s corresponds to lab energies Ti ~ 3 T ~ 60 — 240 MeV. 
Such energies are below the threshold for pion production, and so direct comparison at the 
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relevant energies is not possible. However, one can be very encouraged by the excellent 
agreement between the OPE approximation and the experimental data for Tl = 280 — 1000 
MeV, and the fact that that region is not so far from the energies of interest. Moreover, 
there is no reason to expect a surprise at the lower energies relevant to SN 1987A. Given 
the circumstances, the agreement is much better than one might have expected: OPE does 
not take into account resonances — the threshold for A(1232) production is only Tl — 630 
MeV. Far more detailed treatments of pion bremsstrahlung exist in the literature 15 , in 
which both the A 7 ' (1440) and A(1232) resonances are taken into account; one could improve 
the calculation of axion bremsstrahlung by similar means. However, the energies in the 
supernova are well below the threshold for any baryon resonance, and the additional effort 
seems unjustified. Contrary to the claims of Choi, et al 5 this phenomenological comparison 
seems to validate the OPE approximation. Since the axion emission rate i a oc m 2 , the 
axion mass limit derived scales as eZ 1 ^ — a factor of 2 uncertainty translates into a factor 
of y/2 uncertainty in the mass limit. While it appears that there should be little worry 
about using OPE to compute i a , one must still worry about collective nuclear effects. 
Because of the high densities at the core of the supernova, we have no similar laboratory 
data to compare with. Modulo this important uncertainty, it appears that the rate of axion 
emission from the supernova has been calculated to adequate accuracy, especially given 
the other uncertainties. Likewise, our work seems to justify the use of OPE in calculating 
vv emission from hot neutron stars through the process nucleon-nucleon, neutrino pair 
bremsstrahlung 1 6 . 
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L. Krauss, and F. Wilczek. This work was supported in part by the NSF (at ITP-UCSB), 
the DoE (at Chicago and at Ohio State), and by an Alfred P. Sloan Fellowship. 
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Here we give the exact expression for X^SPIN l-^lx * n the OPE approximation. Recall- 
ing that Mn = (2 m 2 / f n m 2 n )M, it is: 

\M\ 2 = — m 2 ) _2 [( 2 • 4) — m 2 ] 4[(1 • 3) + m 2 ] 

SPIN 

-4mai 3 [(l • a) + (3 • a)] + a 2 3 ^2(1 • a)(3 • a) + m 2 [m 2 - (1 • 3)]) j 
+ above with 1234 — *■ 2143 
4- above with 1234 — + 1243 
-f above with 1234 — > 2134 

-4 (k 2 - - ml)- 1 4m 2 [(l • 2) + (2 • 4) + (2 • 3) - (1 • 4) - (1 • 3) - (3 • 4)] 

+4[— m 4 + (1 • 4)(2 • 3) + (1 • 3)(2 - 4) — (1 - 2)(3 • 4)] 

+2m 3 (ai3 + ai 4 )[(l • a) 4- (3 • a) -J- (4 • a) — (2 • a)] 

+2ma ls ((2 • 3)[— (1 • a) - (4 • a)] + (3 • 4)1(1 • a) + (2 • a)J + (2 • 4)[— (1 • a) - (3 • a)) 

+(1 • 3)[(2 • a) - (4 ■ a)] + (1 ■ 2)[(4 • a) - (3 •<.)] + (1 • 4)((3 ■ a) - (2 ■ «)]) 

+2mi4 ^(2 • 3)[— (1 * a) — (4 * a)] + (3 • 4)[(1 • a) + (2 • a)] + (2 • 4)[ — (1 * a) — (3 • a)] 

+(1 • 3)1(4 • a) - (2 • a)] + (1 ■ 2)[(3 ■ a) - (4 ■ a)] + (1 • 4)[(2 • a) - (3 • a)]) 

+a ls a„ (2m 2 (l ■ a)[(2 • o)-( 3 • a)— (4 • a)]+2(l • a)[-(3 • 4)(2 • a)+(2 • 4)(3 ■ a)+(2 • 3)(4 ■ a)]| 
+m 2 m|[-m 2 + (2 • 3) + (2 • 4) + (1 • 3) + (1 ■ 4) - (3 • 4) - (1 • 2)] 

+m 2 [(1 • 2)(3 • 4) - (1 • 3)(2 ■ 4) - (1 • 4)(2 • 3)]) 

4- above with 1234 — > 1243 
+ above with 1234 — ► 3421 
+ above with 1234 — ► 3412 
+ above with 1234 — » 4312 
+ above with 1234 — ► 4321 
4 above with 1234 — ► 2143 
+ above with 1234 —>2134 
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where we denote the 4- momenta of the incoming and outgoing protons by 1,2 and 3,4 
respectively, the 4-momentum of the pion by a , and k = 2 — 4, / = 2 — 3, a.ij = m/(i • a) + 
m/(j ■ a). 

For consistency with the previous calculations of axion bremsstrahlung we have taken 
the coupling of the exchanged pion to the proton to be 2 m/m*, while, in analogy to the 
axion proton coupling we have taken the coupling of the outgoing pion to be m/ f n . Had 
we taken the neutral pion proton coupling to be either 2 m/m n (or m/f n ) throughout our 
result for cr(pp — ► pp + ir°) would have been larger by a factor of about 2 (smaller by a 
factor of 4). 
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Figure Captions 


FIGURE 1 — The 4 direct and 4 exchange diagrams for nucleon-nucleon, axion (or pion) 
bremsstrahlung, in the OPE approximation. 

FIGURE 2 — A summary of the experimental data for pp — * pp + i r° and the analytical fit 
to the data as a function of the lab KE of the incoming proton (from Refs. 12 and 13). 

FIGURE 3 — The naive OPE prediction for the cross section (using the matrix element 
in Eq.(2), i.e., talcing m n = 0) and the experimental data. In addition, the naive OPE 
prediction as corrected by the approximate threshold correction to the matrix element, cf., 
Eq.(10), is shown. 

FIGURE 4 — The exact OPE prediction for the cross section (using the matrix element 
computed with m K ^ 0), cf., Eq.(ll), and the experimental data. The effect of neglecting 
the 3,3 resonance in the OPE calculation around 600 MeV is apparent. 

FIGURE 5 — The naive OPE cross section (as computed from the matrix element in Eq.(2)) 
and the exact OPE cross section with m K set to zero in the matrix element. Except 
near threshold the agreement between the two is excellent; as described in the text the 
discrepancy at threshold is as expected. 
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